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Abstract. Luminosity and mass functions in the old open cluster NGC 188 are analysed by means of J and H 
2MASS photometry, which provides uniformity and spatial coverage for a proper background subtraction. With 
an age of about 6-8 Gyr, NGC 188 is expected to be suffering the effects of advanced dynamical evolution. Indeed, 
previous works in optical bands have suggested the presence of mass segregation. Within the uncertainties, the 
observed projected radial density profile of NGC 188 departs from the two-parameter King model in two inner 
regions, which reflects the non-virialized dynamical state and possibly, some degree of non-sphericity in the spatial 
shape of this old open cluster. Fits with two and three-parameter King models to the radial distribution of stars 
resulted in a core radius R CO rc = 1.3 ± 0.1 pc and a tidal radius Rtidai = 21 ± 4pc, about twice as large as the 
visual limiting radius. The concentration parameter c = 1.2 ±0.1 of NGC 188 makes this open cluster structurally 
comparable to the loose globular clusters. The present 2MASS analysis resulted in significant slope variations 
with distance in the mass function <j>(m) oc m _ ' 1+x - ) , being flat in the central parts (x = 0.6 ± 0.7) and steep in 
the cluster outskirts (\ — 7.2 ± 0.6). The overall mass function has a slope \ — 1-9 ± 0-7, slightly steeper than 
a standard Salpeter mass function. In this context, NGC 188 is similar to the 3.2 Gyr, dynamically evolved open 
cluster M 67. Solar metallicity Padova isochrone fits to the near-infrared colour-magnitude diagram of NGC 188 
resulted in an age of 7.0± 1.0 Gyr. The best fit, obtained with the 7.1 Gyr isochrone, produced a distance modulus 
(m - M) = 11.1 ± 0.1, E(B - V) = 0.0, and a distance to the Sun d = 1.66 ± 0.08 kpc. The observed stellar 
mass (in the range 0.98 Mq — 1.08 Mq) in NGC 188 is m b s = 380± 12 Mq. A simple extrapolation of the observed 
overall mass function to stars with 0.08 Mq resulted in a total present mass of m to t ~ (1.8 ± 0.7) x 10 4 Mq. On 
the other hand, for a more realistic initial mass function which flattens in the low-mass range, the total mass in 
NGC 188 drops to m to t ~ (3.8 ± 1.6) x 10 3 Mq. Since mass-loss processes such as evaporation and tidal stripping 
have been occurring in this old open cluster for about 7 Gyr, the primordial mass in NGC 188 must have been 
significantly larger than ~4x 10 3 Mq . We also examined the consequences of the presence of unresolved binaries 
and concluded that, even if dominant in the CMD, binaries alone are not responsible for the flat central mass 
function, which supports the mass-segregation scenario. 
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1. Introduction 

With an age of about 7 Gyr (e.g. Carraro & Chiosi ,1994), 
NGC 188 is one of the oldest known open clusters in the 
Galaxy. Located relatively far from the disk, its field is 
neither heavily contaminated by background stars nor ob- 
scured by dust, and contains a few hundreds of member 
stars. These conditions make NGC 188 an excellent labo- 
ratory for testing modern theories of stellar and dynamical 
evolution of star clusters. 
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The main effects related to the dynamical evolution of 
an open cluster are (i) mass segregation, (it) tidal strip- 
ping and disruption by the Galactic gravitational field, 
and (tit) interactions with the disk, since a typical open 
cluster at the solar radius will cross the Galactic plane 
10-20 times before being disrupted (de la Fuente Marcos 
1998). Evaporation of a cluster may result from two-body 
relaxation processes, in which several stars can acquire 
positive energy during the closest approaches and eventu- 
ally leave the cluster. Mass-loss during the course of stellar 
evolution and encounters with giant molecular clouds (e.g. 
Wielen 11991)1 also contribute to reducing a cluster's life- 
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time. Considering all these factors, the destruction time- 
scale for open clusters in the solar neighbourhood is about 
600 Myr (Bergond, Leon, & Guibert 12001(1 . Consequently, 
only the most massive open clusters or those located at 
large Galactic radii are expected to survive longer than a 
few Gyr fFriel 11995(1 . Indeed, de la Fuente Marcos l|1998(l 
has shown that the Milky Way should host several hun- 
dred thousands of open cluster remnants which have been 
disrupted both by the Galaxy and their own dynamical 
evolution. 

Taking into account all of the above factors, as well 
as its location in the Galaxy, the dynamical evolution of 
NGC 188 should be essentially associated to mass segre- 
gation, evaporation and tidal stripping by the Galactic 
gravitational field. 

As a consequence of mass segregation, the central re- 
gions of the more evolved clusters should present a main 
sequence (MS) depleted of low-mass stars, thus creat- 
ing a core rich in compact and giant stars (Takahashi 
& Portegies Zwart GDOOP . Also, the luminosity (or mass) 
function (hereafter LF and MF) of a dynamically evolved 
open cluster should present slope changes with distance 
to the center, being flatter in the central parts than in the 
outer regions. 

Examples of advanced dynamical evolution in star 
clusters have been reported in recent years, e.g. in 
NGC 3680 (Anthony- Twarog et al. lTMTl Bonatto, Bica & 
Pavani l3UM|l and M67 (Bonatto k Bica lMTSjl . in which 
the presence of a core depleted of low-MS stars and a halo 
rich in low-mass stars have been confirmed with 2MASS 
photometry. 

The accurate determination of LF spatial variations, as 
well as the characterization of a cluster's stellar content, 
depends critically on a proper background subtraction. 
The uniform and essentially complete sky coverage pro- 
vided by 2MASS (Two Micron All Sky Survey, Skrutskie 
et al. ITMTjl has proven to be a powerful tool for properly 
taking into account background regions with suitable star 
count statistics (see, e.g. Bonatto Bica & Pavani 2004 
Bonatto & Bica l2"UT)3l Bica, Bonatto & Dutra l2"UU3l Bica, 
Bonatto & Dutra 1277114) . 

In the present work we investigate the dynamical state 
and stellar content of NGC 188. We employ J and H 
2MASS photometry to analyse a large spatial area in the 
direction of the cluster. In particular, we search for spa- 
tial changes in the LF which would characterize the ad- 
vanced dynamical state of this open cluster. This is the 
first analysis of NGC 188 in the near-infrared with wide 
spatial coverage and resolution. 

This paper is organized as follows. In Sect. |2] we provide 
general data on NGC 188 and present a brief review on this 
open cluster. In Sect. we obtain the 2MASS photometry 
and introduce the J x (J — H) colour-magnitude diagrams 
(CMDs). In Sect.0]we discuss the projected radial density 
distribution of stars and derive structural parameters for 
NGC 188. In Sect. [5] we fit isochrones to the near-infrared 
CMD and derive cluster parameters. In Sect. H3 we derive 
the LFs and MFs, estimate the stellar mass still present in 



NGC 188 and examine the presence of unresolved binaries 
in the CMD. Concluding remarks are given in Sect. 
In Appendix E] we discuss the consequences of applying 
colour-filters to the CMD and the background subtraction. 

2. The old open cluster NGC 188 

NGC 188 has been in the spotlight since the early works 
of Sandage (19611 IT962(I when the age of this open clus- 
ter was estimated to be r ago ~ 10 — 16 Gyr, making 
NGC 188 the probably oldest observable open cluster in 
the Galaxy. However, later works revised its age down- 
wards (see Carraro et al. 119941 for a comprehensive age 
census) to the presently accepted value of around 6-7 Gyr 
(von Hippel & Sarajedini [1998 Sarajedini et al. I1999|l . 
The WEBDA open cluster database 1 (Mermilliod 11996(1 
gives T agc w 4.3 Gyr for NGC188. 

The actual, advanced age of NGC 188 still makes it one 
of the oldest surviving open clusters in the Milky Way. 
As such, NGC 188 has been used to test stellar evolution 
models, since its CMD exhibits a wide giant branch (e.g. 
Twarog 1997, Norris & Smith 1985), blue stragglers (e.g. 
Leonard & Linncll 1992, Dinescu et al. 11996(1 and numer- 
ous binaries (von Hippel & Sarajedini 1998, and references 
therein) . Most of the previous research has been restricted 
to optical bands. 

One of the reasons why NGC 188 has survived to such 
an advanced age may be its almost circular, highly in- 
clined external (9.5 < R(kpc) < 11) orbit, which avoids 
the inner disk regions for most of the time ( a plot of the 
Galactic orbit of NGC 188 is shown in Fig. 1 of Carraro & 
Chiosi lT994(l . As a consequence of its orbit, the probability 
of significant dynamical interactions with giant molecular 
clouds becomes considerably reduced and thus, the role of 
internal processes in the dynamical evolution of NGC 188 
increases in importance. 

Based on V and I-band photometry of the central parts 
of NGC 188, (von Hippel & Sarajed ini ITMl , and UBVRI 
CCD photometry, Sarajedini et al. (1099), suggested that 
NGC 188 has a typical initial mass function. But, as a con- 
sequence of the advanced, i.e. mass-segregated dynamical 
state, the low-mass stars presently missing from the cen- 
tral parts are either in the cluster outskirts or have left 
the cluster (von Hippel & Sarajedini 11998(1 . 

WEBDA gives a slightly subsolar metallicity for 
NGC 188, which agrees with other values, e.g. [Fe/H] = 
0.02 ± 0.11 (Caputo et al. HMty . [F e/H] = -0.12 ± 0.16 
(Hobbs, Thorburn & Rodriguez-Bell 1990} and [Fe/H] = 
-0.02 (Twarog, Ashman & Anthonv-Twarog lTMTjl . 

The central coordinates of NGC 188, precessed to 
J2000, are a = 00 h 47™ 53 s and S = +85° 15' 30", which 
convert to £ = 122.85° and b = +22.39°. 

3. The 2MASS photometry 

Despite having lost its title as the oldest Galactic open 
cluster, NGC 188 still deserves attention, since its age im- 

1 http://obswww.unige.ch/webda 
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Fig. 1. J x (J - H) CMDs for the cluster NGC 188 (left 
panel, containing 1804 stars) and corresponding offset field 
(right panel, containing 998 stars). Both CMDs have been 
extracted within equal areas. 



plies advanced stellar and dynamical evolution, and its 
location in the Galaxy favours deep observations, partic- 
ularly in the infrared. We base the present analysis of 
NGC 188 essentially on J and H 2MASS 2 photometry. 
2MASS photometric errors typically attain 0.10 mag at 
J « 16.2 and H « 15.0, see e.g. Soares & Bica (HSU). 
The VizieR 3 tool has been used to extract stars in a cir- 
cular area of 45' radius centered on the coordinates given 
in Sect. |21 In order to maximize the statistical significance 
and representativeness of background star counts, we used 
the stars in the outermost annulus (40' < R < 45') as off- 
set field. See Appendix ^ for a discussion on background 
subtraction. 

In Fig. [2 we show the J x (J - H) CMDs for NGC 188 
(left panel) along with the corresponding (equal area) off- 
set field (right panel). In order to maximize visualization 
of star sequences, the CMD in the direction of NGC 188 
has been built with stars extracted within an area of 15' 
in radius. This dimension is, however, smaller than the 
limiting radius of NGC 188 (Sect.©. 



2 All Sky data release, available 
http:/ / www. ipac. caltech. edu/2mass /releases/ allsky/ 
6 http:/ /vizier. u-strasbg.fr/viz-bin/ VizieR ?-source=II/246 
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The near-infrared CMD of NGC 188 presents a well- 
defined MS with several stars occupying the turnoff, as 
well as a wide giant branch. The large width of the MS, 
which cannot be entirely accounted for by photometric er- 
rors (Sect. EHEIJ, indicates that this cluster may contain 
a large fraction of binaries, confirming previous optical 
observations (von Hippel & Sarajedini 119981 and refer- 
ences therein) . These features are not present in the CMD 
of the comparison field. The nearly vertical sequence at 
0.2 < (J - H) < 0.4 and 10.5 < J < 13.5, present in both 
CMDs but more populated in NGC 188, is probably due 
to disk stars and, to some extent, to the presence of blue 
stragglers in NGC 188 (Leonard & Linnell ITM21 Dinescu 
et al.EHU. 

The number of stars in the CMD in the direction of 
NGC 188 is 1804, while in its comparison field this number 
drops to 998. 

4. Cluster structure 

The overall cluster structure is analysed by means of the 
star density radial distribution, defined as the projected 
number of stars per area in the direction of a cluster, 
shown in Fig.0for NGC 188. 

Before counting stars, we applied a cutoff (J < 16.0) 
to both cluster and offset field to avoid oversampling, i.e. 
spatial variations in the number of faint stars and spurious 
detections. Colour filters ( Appendix |A"|) have also been ap- 
plied to both cluster and corresponding offset field, in or- 
der to account for the contamination by the Galaxy. This 
procedure has been previously applied in the analysis of 
the open cluster M67 (Bonatto & Bica 2003). As a result 
of the filtering process, the number of stars in the direction 
of NGC 188 inside an area of 20' in radius turns out to be 
1201, compared to 583 in its offset field (same area). The 
radial distribution has been determined by counting stars 
inside concentric annuli with a step of 1.0' in radius. The 
background contribution level corresponds to the average 
number of stars included in the external annulus, which 
lies outside the cluster border at 40' < R < 45'. 

The projected radial density profile of NGC 188 is very 
smooth with a well-defined central concentration of stars. 
The relatively large number of stars in this cluster pro- 
duces small Id Poisson error bars along most of the pro- 
file. Considering the profile fluctuations with respect to 
the background level, we can define for NGC 188 a lim- 
iting radius of Rn m = 24 ± 1'. Because of the essentially 
null-contrast between cluster and background star den- 
sity, any statistical analysis for regions beyond Rum would 
produce excessively high Poisson error bars, and, conse- 
quently, meaningless results. Thus, for practical purposes, 
we can consider that most of the cluster's stars are con- 
tained within Rn m . 

First order structural parameters of NGC 188 are de- 
rived by fitting the two-parameter King fTS6Ba| surface 
density profile to the background-subtracted radial dis- 
tribution of stars. The two-parameter King model essen- 
tially describes the central region of normal clusters (King 
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Fig. 2. Projected radial distribution of surface star den- 
sity. The average background level is shown as a narrow 
shaded area; la Poisson errors are also shown. Magnitude 
(J < 16.0) and colour cutoffs have been applied to the 
object and offset field. The dashed line shows a surface 
density profile (two-parameter King model) fit to the ra- 
dial distribution of stars. 



I1966bl Trager, King & Djorgovsky [T555j l. The fit has 
been performed using a non-linear least-squares fit routine 
which uses the error bars as weights. The best-fit solution 
is shown in Fig. as a dashed line, and the resulting core 
radius is R coro = 2.6 ± 0.2'. With a distance to the Sun 
d© = 1.66 ± 0.08 kpc (Sect. |5J|, the linear core radius of 
NGC 188 turns out to be R CO rc = 1-3 ± 0.1 pc and the an- 
gular diameter (2 x Ri im ) of 48 ± 2' converts to a linear 
limiting diameter of 23.2 ± 1.5 pc. 

Within the errors, the observed density profile departs 
from the two-parameter King model in opposite senses, 
particularly for the regions 0.9 < R(') < 3.5 and 3.5 < 
R(') < 8.6 (Fig. EJ). Considering the physical premises of 
the King model, these deviations may be a consequence 
of the non-virialized dynamical state and possibly, some 
degree of non-sphericity in the spatial shape of NGC 188. 

The core radius of NGC 188 is very similar to that of 
the T ag0 = 3.2 Gyr, relaxed open cluster M67 (R CO rc = 
1.2±0.1pc - Bonatto & Bica l200^j) . However, its limiting 
diameter is larger than that of M67 (12.7±0.6pc), proba- 



bly a consequence of NGC 188's older age, more advanced 
dynamical state and larger mass (Sect. l6~2")l . 

4.1. Tidal radius 

The rather massive (Sect. 16 .2f) nature and present large 
galactocentric distance (Sect. [SJ of NGC 188 make it pos- 
sible to fit its projected radial distribution of stars with the 
more complete empirical density law of King (1962 ), which 
includes the tidal radius Rtidai as parameter in the model. 
The tidal radius depends both on the effect of the Galactic 
tidal field on the cluster and the subsequent internal re- 
laxation and dynamical evolution of the cluster (Allen & 
Martos[1988). In particular, Oh & Lin (1992) have shown, 
through numerical simulations, that relaxation affects the 
tidal radius formation. In the case of NGC 188, the tidal 
radius determination is made possible by the spatial cov- 
erage and uniformity of 2MASS photometry, which allows 
one to obtain reliable data on the projected distribution 
of stars for large extensions around clusters. 

The three-parameter model of King l)1962|) is a poten- 
tial source of information both on the internal structure 
of the cluster, by providing reliable determinations of the 
core radius and the central density of stars, and on the 
Galactic tidal field, since the tidal radius is linked to the 
stripping of stars from the cluster by the Galactic tidal 
field. Since the three-parameter King model provides a 
better description of the outer regions of a cluster than 
the two-parameter model, it is particularly suitable for 
deriving a reliable value of Rtidai- 

The best-fit solution, with a correlation coefficient of 
0.96 and an RMS error of 0.83, resulted in R corc = 3.1 ± 
0.2' and Rtidai = 44 ± 9', which convert to R CO rc = 1.5 ± 
0.1 pc and Rtidai = 21 ± 4pc. This solution is shown in 
Fig. |21 for regions more external than R = 8'. The two- 
parameter model fit is also shown in Fig. for comparison 
purposes. 

The present value of Rtidai is nearly twice as large as 
that estimated by Keenan, Innanen & House l|1973J) and, 
within the errors, comparable to the theoretical value cal- 
culated by Allen & Martos l|1988f) . Interestingly, we note 
that the physical structure of NGC 188 extends well be- 
yond the visual limits, since the presently derived value 
of Rtidai is about twice that of Ri; m - We note as well 
that the inclusion of the additional parameter Rtidai in 
the fit produced an increase of about 15% in R CO re, a 
value which differs 2a from that derived with the two- 
parameter King model. However, since the two-parameter 
King model appears to give a better description of the 
inner regions of star clusters (King I1966bl Trager, King 
& Djorgovsky 11995(1 . we will adopt the previous value 
of R C orc = 1-3 ± 0.1 pc as the core radius of NGC 188. 
With the values derived from the three-parameter model, 
the concentration parameter of NGC 188 turns out to be 
c = log f ^Mdai"j = 12 ± 0.1. To check how NGC 188 fits 

V tlcorc / 

in the context of evolved star clusters, we compare its con- 
centration parameter with those of 141 Galactic globular 
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Fig. 3. Fit of King models to the projected radial distri- 
bution of stars of NGC 188. Dashed line: two-parameter 
model; dotted line: three-parameter model. 



clusters (compiled by Harris [1996 ) 4 , which also includes 
the globular cluster sample of Trager, King & Djorgovsky 
(1995}. We conclude that NGC 188 is structurally compa- 
rable to the loose globular clusters, but still more concen- 
trated than « 25% of the sample. 

5. Fundamental parameters 

To maximize cluster membership probability, the anal- 
yses in the following two sections will be restricted to 
stars extracted within Rn m (Sect.QJ. Cluster parameters 
will be derived using solar metallicity Padova isochrones 
from Girardi et al. (20Q2J> computed with the 2MASS J, 
H, and Kg filters 5 . The 2MASS transmission filters pro- 
duced isochrones very similar to the Johnson ones, with 
differences of at most 0.01 in (J — H) (Bonatto, Bica & 
Girardi EE 



Tfj) . The solar metallicity isochrones have been 
selected to be consistent with previous results (e.g. Caputo 
et al-HHOI Hobbs, Thorburn & Rodriguez-Bell UMTl and 
Twarog, Ashman & Anthony- Twarog 1997). For redden- 
ing and absorption transformations we use Ry = 3.2, 
and the relations Aj = 0.276 x Ay and E(J — H) = 
0.33 x E(B — V), according to Dutra, Santiago & Bica 
and references therein. 



The presence of a well-defined turnoff and giant branch 
in the Mj x (J - H) CMD of NGC 188 constrain the 
age solution to the 6.3 Gyr, 7.1 Gyr and 8.0 Gyr Padova 
isochrones. The best fit, shown in Fig. corresponds 
to the T agc = 7.1 Gyr solution. Accordingly, we adopt 
as the age of NGC 188 r ago = 7.0 ± 1.0 Gyr. The Mj 
values are obtained after applying the distance modu- 
lus derived below for NGC 188. From the isochrone fit 
and related uncertainties we derive a distance modulus 
(m - M) = 11.1 ± 0.1, E(B - V) = 0.0 and a distance to 
the Sun d Q = 1.66±0.08kpc. For illustrative purposes, we 
indicate representative stellar masses along the isochrone. 

With the above distance to the Sun, the Galactocentric 
distance of NGC 188 is d GC = 8.9 ± 0.1 kpc. 

At this point it may be useful to estimate the expected 
mass of NGC 188 based on its Galactocentric distance and 
tidal radius. For clusters with a nearly circular orbit we 

'^tidai^ (KingESl, 
Mr, is the 



3M G Q 



i p is the perigalacticon distance and iv± G 



4 http://physwww.physics.mcmaster.ca/ %7Eharris/mwgc.dat 

5 Available at http://pleiadi.pd.astro.it 



can use the relation M c 
where R p 

Galaxy mass inside R p . For the nearly circular orbit of 
NGC 188, R p ks d GC , and at this distance, M G w 1 X 
10 11 Mq (Carraro & Chiosi IT39l| . With these values, the 
expected mass of NGC 188 turns out to be M c ss (3.9 ± 
2.2) x 10 3 M Q . 

6. Luminosity and mass functions 

In this section we analyze the spatial dependence of the 
luminosity and mass functions in NGC 188, which will be 
subsequently used to derive the stellar mass still locked in 
this old open cluster. 

6.1. Luminosity function 

The 2MASS uniform sky coverage has been shown to be an 
ideal tool for studies like the present one, since the entire 
cluster area can be included in the analyses, and represen- 
tative offset fields can be selected around the cluster (e.g. 
Bonatto & Bica l20"U51 Bica, Bonatto & Dutra I2UU3II . As 
a consequence, advanced stages of mass segregation and 
the cumulative effects of the Galactic tidal pull can be 
properly detected and interpreted. 

In Fig. El we show the LFs (0(Mj)) in the J filter 
(shaded area), built as the difference in the number of 
stars in a given magnitude bin between object (continu- 
ous line) and offset field (dotted line). The LFs are given 
in terms of the absolute magnitude Mj , obtained after ap- 
plying the distance modulus derived in Sect.0 The bin in 
magnitude is AMj = 0.25 mag. The LFs in Fig.Elare built 
after applying magnitude (J < 16.0) and colour cutoffs to 
the object and offset field (Sect.[SJ). 

The large projected area of NGC 188, containing ~ 620 
member stars (Sect. 0}, can be used to search for spa- 
tial variations in the stellar content, with statistically 
significant results. Thus, LFs for different regions inside 
NGC 188 have been built following the structures present 
in the radial density profile (Fig.EJ). These LFs are shown 
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Continuous line: star counts in the cluster area; dot- 
ted line: star counts in the offset field; shaded area: 
background-subtracted LF. The turnoff (TO) is indicated 
in each panel. 



6.2. Mass function 



in panels (a) to (d) of Fig. 03 In particular, the LF in 
panel (c) corresponds to the core region, while in panel 
(b) we show the overall (0.0' < R < R lim ) LF. The back- 
ground 0(Mj) has been scaled to match the projected area 
of each region. Representative MS spectral types (adapted 
from Binney & Merrificld 1998) arc shown in panels (c) 
and (d). The turnoff (Mj ss 2.27, spectral type «F3, 
mass « 1.10 M Q ) is indicated in all panels. The excess 
at Mj « 2.5 corresponds to the accumulation of stars 
near the turnoff (Fig. QJ. Hereafter we will refer to the 
maximum number count in the LFs as turnover. 

As expected from the long-term dynamical evolution 
of this old open cluster, the large number of member stars 
contained inside R= 24' (panels (a), (c) and (d)), does 
not present a uniform spatial distribution of MS spectral 
types (magnitude range 3.0 < Mj < 4.5). The region with 
R > 10' seems to have an LF steeper than those of the 
more internal ones, the core region in particular, which is 
characteristic of mass segregation. 



Cluster mass estimates can be made after converting the 
turnoff - turnover LFs into MFs according to 4>(m) = 

0(Mj) 



using the stellar mass-luminosity relation 
taken from the 7.1 Gyr Padova isochrone. To avoid the 
star excess near the turnoff, we restricted this analysis to 
stars with mass < 1.08 Mq (Mj k 2.7). To increase the 
statistical representativeness of our results, we applied this 
method to the core region, 3.0' < R < 10', 10' < R < 24' 
and to the overall LF. The resulting MFs, including la 
Poisson error bars for each data point, are shown in Fig. [HJ 
The spatial variation in the stellar content of NGC 188 
is clearly seen in the MFs shown in Fig. El being flat in the 
core region and very steep in the region 10' < R < 24'. 
We quantify this effect by fitting the function 4>{m) oc 
rn-Ci+x) to the MFs, using a weighted non-linear least- 
squares routine. The results are given in Table ^ The 
number of stars from the high-mass end to just below the 
turnoff (m = 1.08 M Q ) is obtained by numerical integra- 
tion of the LFs in the corresponding mass interval (co. 4). 
Then, multiplying this number by the mass at the turnoff 
gives us an estimate of the mass stored in evolved stars 
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Fig. 6. Mass functions <p(m) restricted to the turnoff- 
turnover region. The MF fits <f>(m) oc m~^ 1+x ' are shown 
as dashed lines. 



( m evoi)' wn i c h is listed in col. 3 of Tabled 111 c °l- 6 we 
provide the MS mass, obtained by numerical integration 
of the MFs from the turnover to m = 1.08 M Q . Col. 7 gives 
the total observed mass contained in each sampled region 
of NGC 188, while the corresponding number of observed 
stars is listed in col. 8. In col. 9 we provide an estimate of 
the observed mass density in each region. Uncertainties in 
the values of x (col. 5) are derived from the least-squares 
fit, and uncertainties in the mass and number of stars come 
from error propagation. 

The advanced, mass-segregated dynamical state of 
NGC 188 is reflected particularly in the flat MF slope in 
the core region (\ — 0.6 ±0.7) compared to the very steep 
slope (x = 7.2 ± 0.6) in the region 10' < R < 24'. Indeed, 
the mass segregation is fully characterized by the positive 
gradient, with respect to the distance to the cluster center, 
presented by the MF slopes (Table QJ. 

At this point it may be interesting to compare the 
dynamical mass of NGC 188 (Sect. with an indepen- 
dent determination of the mass still locked up in stars in 
this cluster by directly extrapolating the derived MF down 
to the theoretical stellar low-mass end mi ow = 0.08 M . 
Thus, according to the MF in Table ^ the total stel- 
lar mass in NGC 188 is m tot - (1.8 ± 0.7) x 10 4 M Q . 
On the other hand, Kroupa, Tout & Gilmore (192} and 



Kroupa (2001) presented evidence that the MFs of most 
globular and open clusters flatten below ~ 0.5 M Q . As a 
consequence, the total stellar mass in NGC 188 is proba- 
bly less than the value derived above, since most of the 
stars are expected to be found in the low-mass range. 
Accordingly, we derive a more conservative total mass 
value for NGC 188 assuming the universal IMF of Kroupa 
l{2TiUT|l . in which x = 0.3 ± 0.5 for 0.08 M Q - 0.50 M , and 
X = 1.3±0.3 for 0.5QM© -0.98 M Q . For 0.98 M - 1.08 M 
we adopt the value derived in this paper, x = 1-93 ± 0.66. 
As expected, the resulting total stellar mass drops to 
mtot ~ (3.8 ± 1.6) x 1O 3 M . This mass value is in close 
agreement with the mass estimated from NGC 188's galac- 
tocentric distance and tidal radius (Sect.EJ. Following the 
procedure described above, we calculate for each region 
of NGC 188 the extrapolated mass, number of stars and 
corresponding mass density. These quantities are listed in 
cols. 10, 11 and 12 respectively of Tabled 

Mass segregation in a star cluster scales with the relax- 



ation time, defined as t r , 



lax 



3 , where t c 



81nN" cross ' 

R/er v is the crossing time, N is the number of stars and 
er v is the velocity dispersion (Binney & Tremaine ^987). 
For NGC 188, R rj 11 pc (Sect.0J, and we assume a typ- 
ical a v « Skins" 1 (Binney & Merrifield HM5|l . Thus, for 
the observed number of stars, t re iax ~ 24Myr, and for the 
extrapolated number of stars, Wax ~ 600 Myr. Both esti- 
mates of t r eiax are much smaller than the age of NGC 188, 
which is consistent with the presence of mass segregation. 
Interestingly, the first estimate of t rc i ax in NGC 188 was 
made by van den Bergh & Sher |1960|) who from star 
counts derived a mass of ~ 900 M Q and a radius of ~ 6.5', 
thus resulting in t ro i ax ~ 64 Myr. 

With an age of about 7Gyr, NGC 188 must have lost 
a significant fraction of its primordial mass through (%) in- 
ternal processes, such as mass loss and disruption through 
the evolution of massive stars in the early phases, and 
evaporation along the cluster's evolution, and (ii) exter- 
nal processes such as tidal shocks in the Galactic disk 
and encounters with giant molecular clouds. During the 
first few million years of a massive cluster's life, as the 
primordial gas is removed from the system, the gravita- 
tional well changes, decreasing the escape velocity. As a 
consequence, the high-velocity stars can be ejected from 
the cluster in a fraction depending on the initial cluster 
configuration, star formation efficiency and velocity distri- 
bution of the stars (Adams 2000, de la Fuente Marcos & de 
la Fuente Marcos 201)2). Subsequently, evaporation sets in 
and starts reducing the cluster mass. Examples of popula- 
tion decline with cluster age, against evaporation, can be 
seen in the numerical models (including stellar mass loss) 
of de la Fuente Marcos & de la Fuente Marcos J2P02). 
Disk shocking plays an important role not only in reduc- 
ing a cluster's mass but also in flattening its spatial shape. 
Indeed, Bergond, Leon, & Guibert f200I| estimate an up- 
per limit of 5-10% for the mass-loss efficiency during a 
single disk crossing, and observe strong flattening in the 
shapes of NGC 2287 and NGC 2548, which they attribute 
to disk shocking. 
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Table 1. Mass- function fit and related parameters. 



Extrapolated 



Region 



(') 
(i) 



(PC) 

(2) 



(M ) 
(3) 



(M ) 
(4) 



(5) 



(M ) 
(6) 



(M ) 
(7) 



(stars) 
(8) 



(M pc->) 
(9) 



(M ) 
(10) 



(stars) 

(ii) 



(M pc- J ) 
(12) 



Core 0.0-1.4 

3.0-10 1.4-4.8 

10-24 4.8-11.6 

Overall 0.0-11.6 



35 
109 

49 
193 



0.95 -> 1.08 
0.98 1.08 
0.95 -t 1.08 
0.98 -> 1.08 



0.6 ± 0.7 
0.8 ± 0.6 
7.2 ± 0.6 
1.9 ± 0.7 



38 ± 3 
105 ± 34 

47 ± 2 
187 ± 12 



73 ± 3 
213 ± 7 

94 ± 3 
380 ± 12 



69 ± 5 
200 ± 5 

91 ± 4 
333 ± 76 



6.7 
0.5 
0.02 
0.1 



609 ± 261 
2260 ± 968 
1090 ± 475 
3760 ± 1610 



2070 ± 1470 
7700 ± 5440 
3760 ± 2670 
12800 ± 9040 



55.5 ± 23.8 
5.2 ± 2.2 
0.2 ± 0.1 
0.6 ± 0.3 



Table Notes. Col. 1: spatial region expressed in arcmin; col. 2: same as Col. 1 in pc; col. 
col. 4: mass range to which the MF has been fitted; col. 5: MF slope; col. 6: turnoff - 



3: mass of stars later than the turnoff; 
turnover mass calculated from the MF 



fit; col. 7: observed mass, sum of m*„ oi and niMs; col. 8: number of observed stars; col. 9: observed mass density; col. 10: 
extrapolated mass; col. 11: extrapolated number of stars; col. 12: extrapolated mass density. 



With respect to NGC 188, the present-day number of 
stars of ~ 10 4 (Table^) is already larger than the primor- 
dial number of stars in the most-populous models available 
in de la Fuente Marcos & de la Fuente Marcos (j2002). 
Consequently, taking into account (i) the sharp popula- 
tion decline with time observed in the populous models of 
de la Fuente Marcos & de la Fuente Marcos lffii02 ). (ii) 
the other mass-reducing processes discussed in the preced- 
ing paragraph and, (Hi) the old age derived for NGC 188 
(~ 7Gyr), it is clear that the primordial population of 
stars in NGC 188 must have significantly exceeded 10 4 
stars. Consequently, its primordial mass must have been 
larger than 4 x 10 3 M Q . 

6.3. Binaries in NGC 188 

The results presented in Sects. 16.11 and 16.21 are based on 
the simple assumption that all stars in the CMD are sin- 
gle objects. However, stars in multiple systems, binaries 
in particular, usually account for a considerable fraction 
of the stellar content in open clusters (e.g. Montgomery, 
Marschall k Janes HMSjl . 

As a consequence of the dynamical evolution in open 
clusters, multiple systems tend to concentrate in the 
central regions, thus changing the initial spatial dis- 
tribution of stars (Takahasi & Portegies Zwart |7000). 
Observationally, the main effect of a significant fraction 
of unresolved binaries in the central parts of a star cluster 
is that the number of low-mass stars is underestimated 
with respect to the higher-mass stars. Thus, the observed 
central LF (or MF) turns out to be flatter than that of 
the actual, resolved stellar distribution. In addition, the 
resulting MF changes will depend both on the binary frac- 
tion and MS mass range (e.g. Kroupa, Tout &; Gilmore 
119911 and references therein). 

With respect to NGC 188, von Hippel & Sarajedini 
(1998) estimate a fraction of « 50% of multiple systems 
in the central region. At this point, the natural question to 
ask is whether this large fraction of binaries in NGC 188 
may account for the observed flat central MF (Fig. H3 and 
Table UJ. 

To address this question we adopted the following ap- 
proach. First we generate a population of single stars ac- 
cording to a standard (x = 1.35) Salpeter MF distribu- 



tion. To be consistent with the 2MASS CMD of NGC 188, 
we used single stars in the mass range 0.15-1. 10 M©. Then, 
according to a pre-defined fraction of CMD binaries, we 
randomly select stars from the original distribution and 
build the specified number of binaries. Mass and lumi- 
nosity biases are avoided in this process (Kroupa [2002 
Mazeh et al. 2005|) and consequently, the probability of 
a star being selected as a binary member depends only 
on the number frequency of its mass. Each star selected 
to become a binary member is subsequently deleted from 
the original population. Based on the 7.1 Gyr isochrone 
mass-luminosity relation, the combined luminosity of each 
pair member is converted back to mass. When all binaries 
are thus formed and their masses stored in the new pop- 
ulation distribution, we add the remaining (not-paired) 
original single stars to the final population distribution. 
The final, observed MF is then built by counting the com- 
bined number of single and paired stars in each mass bin. 
To reach a high level of statistical significance, the initial 
single-star distribution contains 10 6 stars. 

To check how a varying binary population may affect 
the original, single-star MF, we applied the above ap- 
proach to MFs with slopes Ximi 

= 1.0,1.5 and 2.0, and 
for fractions of CMD binaries fbm = 0.5,0.75 and 1.0. 
fbin = 0.75 means that 75% of the observed CMD objects 
are binaries. We illustrate this effect for Ximf = 1.0 in the 
top panel of Fig. \7\ in which the heavy solid line repre- 
sents the MF of the single-star distribution. The final MF 
distributions are shown as dotted, thin solid and dashed 
lines, respectively for the binary fractions fbm = 0.5,0.75 
and 1.0. All MFs in Fig. are scaled so that the total 

number of stars produced by each MF (^J (f>(m) drri^j is 

the same. 

Since most of the single-star distribution is made up 
of low-mass stars, binary formation severely depletes the 
original MF, particularly for m < 0.25 M©. The binaries 
thus formed have a mass ratio q ~ 0.65 and tend to ac- 
cumulate, in varying proportions (according to fbi n and 
Xi m f), in the region 0.3 < m(M©) < 0.9 of the observed 
MF. Thus, considering the distribution of stars with mass 
> 0.25 M©, the final MFs indeed turn out flatter than the 
single-star distributions. However, for stars in the mass 
range from which we derived the present results (mass 
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Fig. 7. Top panel: varying binary fractions are tested for 
a single-star MF of slope XimS = 1-0- The (original) single 
star population MF is represented by the heavy solid line. 
All MFs are scaled so that they produce the same total 
number of stars. Bottom panel: simulated mass-ratio (q = 
M-z/Mx) distribution for MFs with Ximf = 1-0 and 2.0. 
The binary fraction for both distributions is fb; n = 0.75. 

Table 2. Binary effects on MFs 



Ximf 


q 


ibXn 


= 0.50 


fbin 


= 0.75 


fbh, 


= 1.00 


XO.2-1 


1 Xo. 6-1.1 


XO.2-1 


1 Xo. 6-1.1 


XO.2-1 


1 X0. 6-1.1 


1.0 


0.61 ± 0.23 


0.50 


1.12 


0.36 


1.14 


0.26 


1.15 


1.8 


0.65 ± 0.23 


0.90 


1.79 


0.75 


1.88 


0.75 


1.92 


2.0 


0.69 ± 0.22 


1.37 


2.33 


1.25 


2.43 


1.29 


2.50 



Tabic Notes. Ximf * s the original MF slope; xo 2 — 1 
mass range 0.2 < m(MQ) < 1.1 and. xo 6—1 1 
0.6 < m(M ) < 1.1. 



is the MF slope in the 
the slope in the range 



> 0.8 M©), the final distributions become even steeper 
than the original ones. This is particularly true for the 
Ximf = 2.0 MF, since this distribution contains the largest 
fraction of low-mass stars. 

We summarize the above discussion in Table in 
which we give, for each Ximf! the resulting slopes obtained 
from MF fits to the mass ranges 0.2 < m(M ) < 1.1 and 
0.6 < m(M Q ) < 1.1. 

As expected, the 0.2 < m(M ) < 1.1 MF flattening 
degree is largest for the flattest single-star MF, since the 
relative number of higher-mass stars with respect to lower- 
mass stars is favoured in this distribution. Indeed, the av- 



erage binary mass ratio for Ximf — 1-00 is, within the 
uncertainties, the smallest of the MFs tested here. 

Because of the binary formation, the accumulation of 
stars in the observed MFs for mass < 0.9 M Q ends up 
increasing the MF slope in the 0.6 < m(M©) < 1.1 mass 
range. This effect is present for any Ximf and fbin- 

Finally, we present in the bottom panel of Fig. the 
mass-ratio (q = M2/M1) distribution obtained for MFs 
with Ximf = 1-0 and 2.0, derived for stars with masses in 
the range 0.15 — 1.1 M©. The binary fraction in both distri- 
butions shown in Fig. [7] is fbi n = 0.75. Since the mass ratio 
depends only on the masses of the two companions, the 
distributions for different binary fractions are essentially 
the same, but vary with different MF slopes. Interestingly, 
the nearly- flat distribution for Ximf = 1.0 for 0.3 < q < 1.0 
agrees with the results presented in Mazeh et al. ( 2003 ) for 
spectroscopic binaries with primary masses in the range 
0.6 — 0.85 Mq. The similarity between the simulated mass- 
ratio distribution for Ximf = 1-0 (Fig. EJ and that ob- 
served in spectroscopic binaries (Mazeh et al. I2003|l sug- 
gests that Galactic field stars tend to have an MF with a 
slope around x = 1- 

6.3.1. Widening of CMD sequences 

In the magnitude range we are dealing with in this work, 
uncertainties in J and H are typically < 0.10 mag (e.g. 
Soares & Bica |2jQQl2J). Despite the small values, the un- 
certainties produce an additional scatter on the observed 
CMDs and, as a consequence, a widening, particularly of 
the MS. NGC 188 indeed presents a wide MS (Figs. [Hand 
0J, and in this section we investigate how much of the 
MS width can be attributed to binaries and photometric 
errors. The approach here is to produce synthetic CMDs 
which simulate both effects and to compare them with 
actual CMDs of NGC188. 

In 2MASS photometry, the average photometric un- 
certainty (e) for a given magnitude follows a tight relation 
with magnitude, as can be seen in Fig. [HJ in which we plot 
data for NGC 188. For each J magnitude, the photometric 
error is obtained by randomly selecting a value in the in- 
terval (— €j < e < +€j), assuming that e follows a normal 
distribution curve. H errors are obtained similarly. 

We first generate an MS stellar distribution according 
to a given initial mass function (IMF) and binary fraction, 
as described in Sect.EH For NGC 188, the MS is assumed 
to follow the 7.1 Gyr isochrone. For each single-mass star 
and binary thus formed, uncertainties are assigned to the 
J and H magnitudes as described above. Finally, we add 
to the simulated CMD a background stellar distribution. 

We test this approach by comparing the resulting 
CMDs with the observed one for R < 10', shown in panel 
(d) of Fig. [5] the corresponding (same area) background 
field is shown in panel (a). To be consistent with the re- 
sults of Sects. RT^l and l6.3l we use an IMF with x = 0.8 and 
binary fraction fb m = 0.50. In addition, the simulated MSs 
have been generated with approximately the same number 
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Fig. 8. Average photometric uncertainty related to a given 
magnitude. Solid line: en x H; Dashed line: ej x J. Least- 
squares fits result in e H = 0.025 + 1.92 x io-M H /°- 998 ) 
and ej = 0.023 + 4.78 x io-V J / a976 >. 



of stars as the observed, background-subtracted MS. We 
illustrate the individual scatter produced in the IMF dis- 
tribution by the photometric errors and binaries, respec- 
tively in panels (b) and (c). The MSs in panels (b) and 
(c) have similar widths, although, as expected, the pho- 
tometric errors produce scatter symmetrically distributed 
about the isochrone. The binaries, as expected, produce 
scatter biased towards brighter magnitudes. This effect 
was originally pointed out by Kuiper (1935). Finally, as- 
signing photometric errors to the binaries and single stars 
of panel (c), and adding the offset field (panel (a)), we 
obtain the comparison CMDs in panels (e) and (f). For 
illustrative purposes, the 7.1 Gyr isochrone is included in 
all panels. 

The present approach produced near-IR CMD se- 
quence widenings similar to those obtained in the optical 
by, e.g. Kerber et al. ( 2002), Hernandez, Valls-Gabaud & 
Gilmore and Hurley & Tout ijTMgjl . in which bi- 

naries and photometric errors are also included in CMD 
statistical simulations. 

The scatter resulting from the combination of binaries 
(fbin = 0.50) and photometric errors (panel (f)) produces 
a thick MS which, in general terms, presents similarities 
with the observed MS (panel (d)). However, the observed 
MS is clearly wider than the simulated one. Since the pho- 
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Fig. 9. Simulated CMDs compared to the observed one. 
Panel (a) - offset field; (b) - scatter produced by 2MASS 
photometric errors only on a x — 0.8 IMF distribution; 
(c) - scatter produced by a binary fraction fbi n = 0.5 on 
the same IMF of panel (b); (d) - observed CMD in the 
region R < 10'; (e) - CMD resulting from the combination 
of photometric errors, binaries (fbin = 1-0) and the offset 
field; (f) - same as (e) but for fbin = 0.5. The offset field 
distribution is shown by the stars. 



tometric errors are indeed very small in the MS magnitude 
range of NGC 188, the observed excess width is probably 
due to a larger binary fraction. Indeed, a better repre- 
sentation of the observed MS is obtained with fbin = 1-0 
(panel (e)). Interestingly, Kroupa, Tout & Gilmore (1991) 
found that for stars to about 100 pc from the Sun, the LFs 
determined both from photometric parallax and counting 
individual stars are consistent with each other if all stars 
are in binary systems. 

In conclusion, the above arguments show that the pres- 
ence of unresolved binaries in the CMD, even in large 
fractions, cannot account for the large MF flattening de- 
gree detected in the central parts of NGC 188. On the 
other hand, the width of the observed MS can be mostly 
accounted for by a binary population with a fraction 
fbin ~ 1-0, provided that photometric errors are properly 
taken into account. 
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7. Concluding remarks 

In the present work, the physical structure, stellar content 
and dynamical state of the old open cluster NGC 188 have 
been analysed. The spatial dependences of the luminosity 
and mass functions were studied in detail. The present 
analyses made use of J and H 2MASS All Sky data release 
photometry. 

NGC 188 is about 7 Gyr old and has survived past the 
age at which most open clusters have already dissolved, 
possibly leaving behind a cluster remnant (Carraro 2002 
Pavani et al. 2002 2003). One of the reasons for such a 
long life may be the orbit of NGC 188, which avoids the 
inner disk regions most of the time, thus decreasing the 
probability of large-scale encounters with giant molecu- 
lar clouds. The large Galactocentric distance of NGC 188, 
dec ~ 8.9 kpc, also accounts for the long-time survival of 
this open cluster, together with its relatively large mass. 

NGC 188 presents a smooth radial distribution of stars 
(Sect. 0|, with small la Poisson error bars because of the 
large number of member stars. Its projected radial den- 
sity profile has a marked central concentration of stars 
(Fig- 13 ■ From a two-parameter King model fit to the pro- 
jected radial density profile we derive a core radius R CO rc = 
1.3 ± 0.1 pc. The tidal radius, Rtidai = 21 ± 4pc, has been 
derived by fitting a three-parameter King model to the 
profile. The tidal radius is nearly twice as large as the visu- 
ally determined linear limiting radius of 11.6±0.8pc. The 
concentration parameter of NGC 188, c = 1.2±0.1, makes 
this open cluster structurally comparable to the loose 
globular clusters, but still more concentrated than « 25% 
of them. The observed projected radial density profile of 
NGC 188 departs from the two-parameter King model in 
two inner regions, which reflects the non-virialized dy- 
namical state (despite the cluster's old age) and possibly, 
some degree of non-sphericity in the spatial shape of this 
old open cluster. Its Mj x (J - H) CMD (Fig. EJ, with 
a well-defined turnoff and giant branch, as well as evi- 
dence of blue stragglers and binaries, can be fitted with 
the 7.0 ± 1.0 Gyr solar metallicity Padova isochrones. The 
best-fit with the 7.1 Gyr isochrone and related uncertain- 
ties result in (m - M) = 11.1 ± 0.1, E(B - V) = 0.0 and 
d Q = 1.66 ±0.08 kpc. 

Mass segregation in NGC 188 is reflected in the spatial 
variation of its MF (Fig.0), which is flat (\ = 0.6 ±0.7) in 
the core region and steep (\ = 7.2 ± 0.6) in the outskirts. 
The MF fit <j)(rn) oc m~( 1+x ) resulted in an observed stel- 
lar mass (MS and giants) of « 380±12 M . If one extrapo- 
lates the 2MASS MF fit down to the theoretical low-mass 
end mi ow = 0.08 M© the total stellar mass in NGC 188 
turns out to be ~ (1.8 ±0.7) x 1O 4 M . However, if we as- 
sume a more representative IMF, which flattens for masses 
below ~ 0.5 Mq, the total stellar mass in NGC 188 turns 
out to be ~ (3.8 ± 1.6) x 1O 3 M . This value is in close 
agreement with the mass estimated for a tidally truncated 
cluster at the present position and with the same tidal ra- 
dius as NGC 188. The primordial mass in NGC 188 must 
have been significantly larger than ~ 4 x 10 3 M Q , since 



mass-loss processes such as evaporation and tidal strip- 
ping have been acting on this cluster for about 7 Gyr. 

We also addressed the observational consequences of 
the presence of unresolved binaries in the central parts of 
NGC 188, since these systems may dominate the 2MASS 
CMD. After testing the effects of different binary fractions 
and single-star distributions, we conclude that unresolved 
binaries alone cannot account for the flatness of the cen- 
tral MF in NGC 188, thus supporting the mass-segregation 
scenario. Finally, the large width exhibited by the MS of 
NGC 188 cannot be explained only by photometric errors; 
instead, it is consistent with a binary fraction of nearly 
100%. 
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Appendix A: Effects of colour filters and 
background selection 

Since star clusters are usually projected against relatively 
rich Galactic fields, the observed CMDs may be heavily 
contaminated by background stars, which thus introduce 
additional noise in the photometry. This effect is partic- 
ularly critical for objects projected against the Galactic 
disk or bulge, fn this context, a careful subtraction of 
this contamination, usually by means of colour and mag- 
nitude filters, is necessary. In this section we will examine 
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Fig. A.l. The dashed line shows the colour filter used for 
NGC 188. Magnitude cutoffs are J = 5 and J = 16. 



the consequences of having applied filters to the CMD of 
NGC 188, both in its physical structure and LF. 

NGC 188 lies at b = +22.39°, consequently its CMD 
is not heavily contaminated, except for disk stars (nearly 
vertical sequence at (J — H) w 0.2 — 0.4) and the faint ones 
for J > 16. Thus, we used a colour filter which follows the 
7.1 Gyr isochrone shape and preserves most of the stars. 
The filter is shown as a dashed line in Fig. IA.1I 

In Fig. IA.2I (top panel) we compare the radial distri- 
bution of stars obtained with no filter (circles) and that 
using the colour filter (squares). The background level 
has been subtracted from both curves. The filtering ef- 
fect, besides decreasing the number of stars, produces a 
thinner and smoother profile with respect to that without 
filter. Indeed, R CO rc = 2.62' for the filtered profile while 
R-corc = 3.14' for the unfiltered one. 

The overall (0' — 22') filtered and unfiltered LFs are 
shown in the bottom panel of Fig. IA.2I The two curves 
are very similar from J — 8 to J ss 13.5. As expected, 
the differences begin to increase for fainter magnitudes, in 
particular for J > 14.5, the magnitude range where most of 
the discarded stars are found. If these stars were preserved 
in the analysis, they would produce an abnormally steep 
profile which would artificially increase the number of low- 
mass stars in the cluster. 
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Fig. A. 2. Top panel: background-subtracted radial distri- 
bution of stars. Bottom panel: overall LFs. 



We conclude that the use of a colour filter on a CMD 
does not introduce significant bias or artifacts either in the 
radial distribution of stars or in the LF of a star cluster, 
at least for those clusters resembling NGC 188. Besides, 
the changes introduced by filters are expected to occur in 
the right direction, since they account for the background 
contamination. 

Background selection and subtraction is another criti- 
cal step in the analysis of star clusters. Ideally, the offset 
field should represent the background stellar contribution 
to the cluster field, both in number and spectral type dis- 
tribution. This is crucial in particular for highly dynami- 
cally evolved clusters in which a significant fraction of the 
low-mass stars are dispersed in the external regions. For 
these clusters the offset field selection should be guided 
by a compromise between distance from the cluster bor- 
der and representativeness of the background stars. 

The uniform sky coverage of 2MASS provides a means 
to suitably account for the stellar background. We illus- 
trate this process in Fig. IA.3I in which we show the cu- 
mulative distribution of stars for three circular, external 
regions of NGC 188, 30' - 35', 35' - 40', 40' - 45'. The first 
two distributions have been scaled to match the area of 
the outermost distribution. 

Within the number fluctuation represented by the 
Poisson error bars, the three distributions are similar, 
which reflects the near uniformity of the stellar content 
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Fig. A. 3. Cumulative distribution of stars in external re- 
gions of NGC 188. Dashed line: 30' - 35'; Dotted line: 
35' — 40'; Solid line: 40' — 45'. Poisson error bars for the 
outermost distribution are shown. 

around NGC 188. However, more internal distributions 
present a slight excess of low-mass stars with respect to 
the outermost distribution, which may be accounted for 
by the long-term effects of mass segregation near the tidal 
radius. For this reason we selected the region 40' — 45' to 
represent the stellar background for NGC 188. 



